Acoustic Emission from Paper Fracture 
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We report tensile failure experiments on paper sheets. The acoustic emission energy and the 
waiting times between acoustic events follow power-law distributions. This remains true while the 
strain rate is varied by more than two orders of magnitude. The energy statistics has the exponent 
f3 ~ 1.25 ± 0.10 and the waiting times the exponent r ~ 1.0 ± 0.1, in particular for the energy 
roughly independent of the strain rate. These results do not compare well with fracture models, for 
(brittle) disordered media, which as such exhibit criticality. One reason may be residual stresses, 
neglected in most theories. 
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A simple, common-day statistical physics experiment 
is to tear a piece of paper into two parts. During this 
noise is heard []], evidently originating from damage to 
the paper or the propagation of a crack. It is easy to 
make the fracture surfaces "rough" as told by the naked 
eye without any more sophisticated analysis, and the ex- 
perience is that it is difficult to do the tearing so that 
the cut is at all close to straight. Meanwhile, one can 
confirm the existence of disorder in the "sample" by just 
looking through the sheet against any reasonable light 
source. The paper looks cloudy, which comes to a large 
degree from the local fluctuations of areal mass density. 
This tells that elastic and fracture properties vary locally. 
Such a test poses the related questions: how does a crack 
develop in an inhomogeneous material, and what kind of 
properties does the acoustic emission have? 

In fracture there is evidence of scaling properties fa- 
miliar from statistical physics. The roughness of cracks, 
as measured by e.g. the root mean square height fluctu- 
ations vs. length scale, can often be described by self- 
afline fractal scaling [3 C3 Q Another such quantity is 
the distribution of strength, and its average versus sam- 
ple size. In particular for (quasi)-brittle materials, the 
average strength can follow scaling laws, if the sample 
size is varied, which derive from the presence of disorder 
in the material HIE0- 

Investigating such questions combines materials sci- 
ence and statistical mechanics. Two fundamental prob- 
lems are, what is the role of disorder in fracture? And, 
how do cracks develop and interact in real materials? 
These join forces in the failure of a notched three- 
dimensional sample, visualized by the advancement of 
an one-dimensional crack front or line in the presence of 
disorder 0, 0. It has an equation of motion, in an en- 
vironment with varying properties as local strength and 
elastic modulus. The trail left by the line- like front is the 
two-dimensional crack surface, which can become rough 
(self-afBne) if the propagating line-like object develops 
critical fluctuations. However, it is difficult to propose 
such a model that would match with the experimental 
data. For one, simple models tell that the disorder and 



microcracks interact in crack dynamics, so that some- 
times it is not correct to consider an isolated object like 
the front or a crack tip. From the materials research 
viewpoint the dynamics is interesting since it deals with 
the engineering quantities of strength and toughness, or 
the resistance to intrinsic flaws. For instance in fiber and 
two-phase ceramic composites one can control these two 
by varying the mixing and constituents. 

We study fracture in ordinary paper, as a two- 
dimensional disordered material, via acoustic emission 
(AE) analysis. The release of acoustic energy is related in 
paper to irreversible deformation, microcracks, and, per- 
haps, to plasticity. A large-scale analogy is earthquakes. 
Their energies are described by a power-l aw p robabil- 
ity distribution, the Gutenberg-Richter-law |l(J with an 
energy exponent (3. Cracks in paper may be self-affinc 
[111 E2], with a roughness exponent close to that in a 
scalar two-dimensional fracture model, random fuse net- 
works (RFN's) [H E3- Both these are near 2/3, the 
one for surfaces of minimal energy f [it, IeH llij). Paper 
provides an example of crack growth in the presence of 
disorder. The probability distributions of the released 
energy and the temporal statistics, describing dynamical 
aspects, can be compared to models for two-dimensional 
failure, and to other work on AE. 

Here the strain rate e is varied in the tensile test, to 
study possibly time-dependent effects in the fracture pro- 
cess. Our first main conclusion is that in strain-controlled 
tests power-laws are found in the statistics, with cither 
no or at most a weak dependence on the strain rate. The 
power-laws do not follow the predictions of theoretical 
models outlined below, as fuse networks or mean-field 
ones. It is als o ap parent that an eventual localization of 
the crack 0, 0] does not change the microscopic prop- 
erties of the crackling noise, as long as the crack prop- 
agation is stable. This implies that paper failure is not 
a "phase transition" with a diverging correlation length, 
but leaves open the origins of the observed scalings. 

To describe the interactions in an elastic medium 
with randomness incorporated one starts from mean-field 
stress-sharing. In fiber bundle models (FBM) the applied 
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external force F(t) is shared by all the N(t) fibers, and as 
the (random) failure threshold of the currently weakest 
one is reached, the stress a(t) = F(t)/N(t) is instanta- 
neously distributed among the remaining fibers. As a 
sign of criticality there is a divergent scale, the size of 
a typical avalanche which is made up of all the fibers 
that break down due to a single, slow-time scale stress- 
increment. The total avalanche size (N) distribution fol- 
lows P(N) ~ N~ 5 / 2 [HE3- Thus S lobal load sharing 
fracture resembles a second-order phase transition. 

Local stress enhancements can be added to FBM's 
e.g. by considering fiber chains in which the stress from 
microcracks, of adjacent missing fibers, is transfered to 
the ones neighboring the crack. A more catastrophic 
growth results (with an exponent much larger than 5/2), 
resembling a first-order phase transition since the elas- 
tic modulus has a finite drop at a c /e c 0, 0, |2£j. A 
finite-dimensional, but scalar, approximation is given 
by random fuse networks [2]]]. RFN's reproduce with 
strong enough disorder many of the features of mean- 
field FBM's, in spite of including the stress-enhancements 
and shielding from other microcracks |U HH . The dy- 
namics of crack growth is not understood well even in 
this simple model, including why the cracks seem to 
become self-affine, with a roughness exponent close to 
2/3 14]. The models' event statistics can be compared 
with experiments by considering the energies. For the 
RFN, in an event the lost energy is Ei ~ AGief, where 
AGi ~ Nbroken is the change in the "elastic modulus" 
due to the AE event, at e,. Simulations on RFN's for 
strong disorder that reproduces the FBM exponent (5/2) 
for the avalanche sizes yield the exponent [3rfn ~ 1.8 
and the FBM value is Pfbm ~ 2 0. 

Experiments on, mostly, three-dimensional systems 
[H El IHIH H3 have yielded power-laws for the AE 
energy release; the typical exponent (3 is 1 . . . 1.5. One 
idea is that approaching final failure resembles a phase 
transition: in the Lyon group's stress-controlled exper- 
iments indications were seen of a critical energy release 
rate 0] . This means that the AE energy would diverge 
like a power-law in the proximity of a critical point, the 
sample strength a c . Another problem is the nature of 
the disorder in the material: whether quenched disorder 
is able to give such power-law statistics for the energy 
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Normal newsprint paper samples were tested in the 
machine direction [2^ on a mode I laboratory testing ma- 
chine of type MTS 400/M. Due to the lack of constraints 
the samples could have out-of-plane deformations, and 
none of the three fracture modes (I, II, III) is excluded 
on the microscopic level. The deformation rates e var- 
ied between 0.1 %/min and 100 %/min. The AE system 
consisted of a piezoelectric transducer, a rectifying ampli- 
fier and continuous data-acquisition. The time-resolution 
of the measurements was 10 fis and the data-acquisition 
was free of deadtime. The stress was measured simulta- 



neously to AE with a time resolution of 0.01 s. We made 
20 identical repetitions for statistics. The 100 by 100 mm 
samples had initial notches (size 15 mm) to achieve stable 
crack growth. Typical sheet thickness is about 100 /im. 
The fracture statistics are not affected by such a notch 
lii). For these strain rates the sound velocity timescale 
is much faster than that implied by 1/e. Each individual 
test contributes, at most, 1000 - 2000 events, so we can 
only look at integrated probability distributions and not 
in detail at local averages of e.g the acoustic event size 
vs. e. 

Fig. JQ| shows an example of two tests under strain- 
control. Stress-strain curves have typically three parts: 
pre-failure (almost) linearly elastic one, the regime close 
to the maximum stress where the the final crack starts 
to propagate or is formed, and a tail that arises due to 
the cohesive properties of paper which allow for stable 
crack propagation. The faster the strain rate the less is 
the role of the tail. For the smallest strain rate most 
of the AE originates from tail (more than 90 %) while 
for the highest the situation is the opposite. Quantities 
of interest are the statistical properties before the max- 
imum, after it, and the integrated totals, in particular 
the energy distributions. The time series of events allows 
to make qualitative observations of the correlations be- 
tween subsequent events and to draw conclusions about 
the event properties as such. For this rather brittle pa- 
per grade, and the strain rates used, the elastic modulus 
is independent of e, and by AE we are able to detect a 
constant fraction of the elastic energy j24j . 

Fig. J5J shows the scaling of the energy for a fixed 
strain rate. The behavior is power-law -like, with sev- 
eral orders of magnitude of scaling. The same exponent 
fits all three different cases: the pre and post maximum 
stress cases, and the sum distribution. If the strain rate is 
varied the same conclusion holds and the exponent only 
fluctuates at random [24|- We have j3 = 1.25±0.10, in dis- 
agreement with the fiber bundle ones and with that from 
the fuse networks, though 0rfn is closer than /3fbm- 
The practical implication is, that the material can with- 
stand more damage than expected since j3 < Pmodeis- In 
paper, the post-maximum events correlate with the ad- 
vancement of the final crack. The remaining ligament 
length/width contracts, thus the elastic modulus drops 
should (assuming a constant stress state) here be quali- 
tatively related to AE event energies (see also US)- 

In an experiment with a varying strain rate e both 
the event durations and the waiting time between any 
two events may (consider Fig. Q) depend on e. If the 
crack dynamics becomes "fast" then events take place on 
timescales set by the sound velocity. This establishes a 
timescale t s — Ax/v s where Aa; is the spatial separation 
of two events and the sound velocity » s ~2x 10 3 m/s. In 
a sample of linear size 0.1 m this results in a maximum 
t s ~ 10~ 4 s. In the failure of brittle carbon foams the 
eventual critical crack growth is dictated by such fast 
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events [lij . 

In quasi-static fracture models the dynamics of cracks 
is assumed such that the stress field is equilibrated in- 
finitely fast during microfracture events, between further 
adiabatic increases of strain or stress. Thus the only 
time-dependence of any temporal statistics is in the av- 
erage time interval, between AE events, proportional to 
1/e. For the FBM model there are no correlations be- 
tween the waiting times and event size or durations, ex- 
cept for the trend that the average waiting time decreases 
as <r c is approached. It is 'critical' and thus the energy 
release rate follows a power-law close to a c . Waiting time 
results do not generally exist for more complicated mod- 
els ■ One can compare to the experimental signatures 
of the intervals between AE events and the integrated 
energy release rate integral J deE event . As noted above 
there is some evidence - from stress-controlled experi- 
ments - for a critical behavior for this quantity |l5|. 

Fig. Q demonstrates waiting time r distributions for 
different strain rates. There is a clear power-law, whose 
exponent (t) remains roughly the same for all the strain 
rates [3l|. Importantly, this is true for the post-events 
regardless of the origin (before/after <r c ) of the major- 
ity of the AE energy. For the pre-events there might be 
some evidence of the exponent increasing with e. In the 
time-series of events, those with long durations are sepa- 
rated, on the average by shorter intervals from the neigh- 
boring events before/after. Fig. @ depicts the waiting 
times prior to an event with two different data analysis 
methods for distinguishing between possibly overlapping 
events [i^. The interval separations are similar for both 
the post/pre-phases, implying that the microscopic fail- 
ure dynamics does not differ, and that a c or e c can not be 
inferred from the event characteristics. It may be so that 
events which are relatively long are precluded by longer 
waiting times. The durations of the events St and the 
sizes are roughly power-law related as (E) ~ (St) 3 . 

The integral J deE event demonstrates a rapid exponen- 
tial growth above a typical strain of e ~ 0.5 %. This orig- 
inates from increasing event sizes, not from an increasing 
density vs. strain. For e > 0.6 % the samples start to fail. 
In the regime where the exponential growth takes place 
the samples develop plastic, irreversible strain with 
a roughly exponential dependence on strain |2 l| . This 
does not imply that the AE measures plastic deformation 
work, mostly, since the rate of increase of e p i is much less 
than that of the energy integral. The exponential strain- 
dependence of the AE implies a typical lengthscale (as 
should be true for the development of plasticity), which 
in turn should be related to crack localization. 

Concluding, failure of ordinary paper shows several 
features associated with critical phenomena. Our take 
on the experiment is that it show that i) there is no clear 
sign of a "critical point" , or a phase transition, in spite 
of the fact that the material is close to linearly elastic. In 
particular, while the event intervals and the event ener- 



gies follow power-law -like statistics, not all the quantities 
do so. Also, it) the temporal behavior hints of compli- 
cated time-dependent phenomena not directly related to 
the fast relaxation of stress. A possible candidate is the 
viscoelastic nature of the wood fibers in paper, but note 
that the macroscopic stress-strain behavior remains al- 
most linearly elastic while AE events already occur. Wc 
suspect the gradual release of internal stresses plays a 
role, perhaps due to frictional pull-out of fibers from the 
network. Finally, Hi) the power-laws as obtained are off 
those predicted by simple fracture models. It remains 
to be seen whether these models can be tailored closer 
to such tensile experiments |33j |. One suggestion is that 
the dynamics of energy release during the events follows a 
different course from the model rules, in particular finite- 
rate dynamics allows for stress overshoots [3 El • The 
relation of the acoustic emission to why cracks get rough 
may be indirect. The latter could relate to the develop- 
ment of pre-failure plastic deformation. 
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FIG. 2: The three (post/pre/total) energy distributions, e 1.0 
[%/min] (for which the post/pre -contributions to AE energy 
are roughly the same). 




FIG. 3: Plots of the event interval r (waiting time) distribu- 
tions, for three strain rates. 
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FIG. 4: The waiting time as a function of event length At, e 
1.0 [%/min]. 
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FIG. 5: Increase of the integrated AE energy vs. strain e. 



